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The microreversibility principle implies that the conductance of a two-terminal Aharonov-Bohm 
interferometer is an even function of the applied magnetic flux. Away from linear response, how- 
ever, this symmetry is not fulfilled and the conductance phase of the interferometer when a quantum 
dot is inserted in one of its arms can be a continuous function of the bias voltage. Such magne- 
toasymmetries have been investigated in related mesoscopic systems and arise as a consequence 
of the asymetric response of the internal potential of the conductor out of equilibrium. Here we 
discuss magnetoasymmetries in quantum-dot Aharonov-Bohm interferometers when strong electron- 
electron interactions are taken into account beyond the mean-field approach. We find that at very 
low temperatures the asymmetric element of the differential conductance shows an abrupt change 
for voltages around the Fermi level. At higher temperatures we recover a smooth variation of the 
magnetoasymmetry as a function of the bias. We illustrate our results with the aid of the electron 
occupation at the dot, demonstrating that its nonequilibrium component is an asymmetric function 
of the flux even to lowest order in voltage. We also calculate the magnetoasymmetry of the current- 
current correlations (the noise) and find that it is given, to a good extent, by the magnetoasymmetry 
of the weakly nonlinear conductance term. Therefore, both magnetoasymmetries (noise and con- 
ductance) are related to each other via a higher-order fluctuation-dissipation relation. This result 
appears to be true even in the low temperature regime, where Kondo physics and many-body effects 
dominate the transport properties. 

PACS numbers: 73.23.-b, 73.50.Fq, 73.63. Kv 



I. INTRODUCTION 

Transport in electric conductors is governed by funda- 
mental principles when the fields applied to the system 
are small. For instance, in the linear regime microscopic 
reversibility leads to symmetric response coefficients, as 
demonstrated by Onsageri In the case of a conductor 
coupled to two terminals, the linear conductance Go is 
an even function of the magnetic field B£ When the con- 
ductor is reduced to typical sizes less than the phase- 
breaking length, electron mesoscopic transport depends 
on the particular arrangement of the attached probes in 
a multiterminal configuration in such a way that current 
and voltage terminals must be exchanged to recover the 
Onsager symmetry^ On the other hand, in the meso- 
scopic regime interference effects associated to the wave 
nature of carriers can be detected in a transport mea- 
surement. A prominent example is an Aharonov-Bohm 
interferometer with a quantum dot inserted in one of its 
arms for which the linear conductance Go is periodically 
modulated by the externally applied flux. However, the 
conductance phase <5, which can be related to the trans- 
mission phase through the quantum dot, shows abrupt 
jumps as a function of the gate voltage^ since the On- 
sager symmetry establishes that S can be or it only- 5 
Further theoretical^ and experimental works&2 have ad- 
dressed the effect of electron-electron interactions inside 
the quantum dot. 

Away from linear response, the principle of microscopic 
reversibility is, generally, not satisfied and, as a conse- 
quence, the two-terminal current, which consists of lin- 



ear as well as nonlinear coefficients, is not a symmetric 
function of B. Recently, this magnetoasymmetry effect 
has been theoretically demonstrated-^— and experimen- 
tally verified22r— . Magnetoasymmetries arise because 
the charge response of the system is, generally, not sym- 
metric when the field orientation is invertedi 1 ^ 11 Out of 
equilibrium, the piled-up charge injected from the exter- 
nal reservoirs is partly balanced by the screening poten- 
tial of the conductor. This internal potential is not an 
even function of B (the Hall potential is a paradigmatic 
example)^ -, leading to magnetoasymmetries seen already 
in the secod-order coefficients within an expansion of cur- 
rents in powers of voltages. 10 Now, computation of the 
internal potential due to long range Coulomb interaction 
is a difficult task and requires self-consistency. This cal- 
culation can be achieved within a mean-field scheme, as 
previous works have done^°-~— However, the importance 
of strong electron-electron correlations such as those giv- 
ing rise to the Kondo effec t 31 > 32 has not been clarified 
yet. This is the goal we want to accomplish in this work. 

Our calculations are also relevant in view of recent de- 
velopments that relate the magnetoasymmetries of the 
current and that of the noise to leading order in a voltage 
expansion^ 3 - - — It has been shown that novel fluctuations 
relations hold in the weakly nonlinear regime between the 
asymmetric second-order conductance and the first-order 
noise susceptibility in terms of a higher-order fluctuation- 
dissipation theorem. These works explicitly check this 
relation for specific systems by treating interactions in 
the mean-field limit— ~— Thus, it is highly desirable to 
find systems in which the nonequilibrium fluctuation re- 
lations can be checked beyond the mean-field case. 
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FIG. 1. (Color online) Sketch of a mesoscopic interferometer 
with a quantum dot inserted in the lower arm. 

We here consider a quantum dot embedded in one of 
the arms of a two-terminal mesoscopic interferometer. 
We employ the nonequilibrium Keldysh Green function 
formalism to describe the transport properties of the sys- 
tem and use an equation-of-motion technique to calculate 
both the current and the noise in the nonlinear regime 
including electron-electron correlation terms not present 
in the Hartree-Fock (mean-field) approximation. 

The oscillations in the conductance of an Aharonov- 
Bohm ring occur in the mesoscopic regime, where the 
electrons' phase coherence is preserved. Interference 
takes place between electron waves that pick up differ- 
ent phases while traversing the two arms of the interfer- 
ometer even if the electron does not directly experience 
the flux enclosed by the ring. When a finite bias is ap- 
plied, the phase of the differential conductance shows a 
continuous variation with B since the Onsager symme- 
try need not hold away from linear response ^ 39 ' Such 
phase rigidity breakings have been recently related to 
the onset of inelastic cotunneling of a Coulomb-blockaded 
quantum dot placed in one of the armsj 41 ' 42 The effect is 
rather generic as the phase symmetry can be also broken 
using microwave fields^ 3 - or coupling to phonons^ 4 - 

We now give a simple argument that sheds light on 
the appearance of magnetoasymmetries in a quantum- 
dot Aharonov-Bohm ring. A sketch of the system is 
depicted in Fig. [TJ We model the interaction in the 
dot with an on-site charging energy U. We assume for 
simplicity that the dot contains a single energy level Sd 
which acquires a finite lifetime 2H/T due to coupling to 
the external reservoirs. Then, within the Hartree approx- 
imation for the Anderson Hamiltonian 4 ^ the (retarded) 
Green function for dot electrons with spin a reads, 



w - e d - u{n ds ) + iVTTrcos (f ) + §r ' 

(1) 

where T renormalizes T due to the presence of the non- 
resonant channel (the upper arm in Fig. [1} . Its (energy- 
independent) transmission is denoted with 71- Impor- 
tantly, the renormalization of the level £d is given by two 
terms. The term that explicitly depends on 71 is an even 



function of the magnetic flux $ and it fulfills the Onsager 
symmetry. However, the term U{nds), proportional to 
the interaction strength, is generally not symmetric un- 
der field reversal since the dot occupation {rid&) out of 
equilibrium need not fulfill the Onsager symmetry. 

Let us expand the dependence of {rids) on the external 
bias V in powers of V, 

(n dff ) = ( nd9 )W + (n d9 ) W V + 0{Vf . (2) 

{nda)^ is the equilibrium charge and must be B- 
symmetric. The nonequilibrium response of the dot to 
leading order in V is given by {rids)^ = d(nda) /dV\v=o- 
{n>do)^ is then determined from the change in the dot 
occupation when a small shift is applied to the leads' 
electrochemical potential. Hence, (nda)^' is a charge 
susceptibility that includes information about the screen- 
ing properties of the dot and, as such, must be com- 
puted self-consistently in the presence of V (e.g., from 
charge-neutrality condition).— We below show that pro- 
cesses of charge filling of the dot from the left lead con- 
tribute to the occupation with a term proportional to 
l+VVfe sin(e<I>/?j) whereas the contribution from the right 
lead is proportional to 1 — \/% sin(e$/fi,). Clearly, the 
sine terms are not even under B- reversal. As a result, 
electron transfer from left to right at a given orienta- 
tion of B does not occur with the same probability that 
the reverse transfer when the B direction is inverted. In 
other words, the injectivity from the left, which is the 
partial density of states associated to carriers injected 
from the left contact^ does not equal the right injec- 
tivity under B reversal* 1 ^ As a consequence, one finds 
(n dg ) {1) oc v / 76sin(e$/fi). Inserting Eq. © in Eq. ([TJ), 
there arises in the denominator a term proportional to 
sin(e$//i), which is responsible, to leading order, 
for the magnetoasymmetry of the nonlinear conductance. 
Note that this term vanishes in the absence of interac- 
tions (U — 0) or at equilibrium (V — 0) in which cases 
the Onsager symmetry is recovered. This further demon- 
strates that magnetoasymmetries arise as a consequence 
of the presence of both interactions and external bias. 

The asymmetric behavior discussed here can be veri- 
fied experimentally. Recent experiments with rings show 
microreversibility violations in the nonlinear regime. An 
unexpected even-odd behavior has been revealed by 
Leturcq et al^ In an expansion of the observed current- 
voltage characteristics, they find that the odd (even) co- 
efficients are symmetric (asymmetric) under reversal of 
B. This is surprising since one would expect all coeffi- 
cients beyond the linear response to be asymmetric, not 
only the even ones. Therefore, it is instructive to derive 
the conductance series expansion. The magnetoasym- 
metric effect appears in the nonlinear regime only be- 
cause the internal potential is an asymmetric function of 
B at finite bias. And this dependence on the internal 
potential is shown in the even coefficients, not in the odd 
ones. A mean-field description gives this behavior— since 
the dot potential depends linearly with V, as can be also 
inferred from our discussion above. Below, we find that 
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the effect persists even if higher electronic correlations 
are taken into account. 

Microreversibility at linear response also leads to the 
fluctuation-dissipation theorem, which relates the dissi- 
pative part of the electric transport (the linear conduc- 
tance) to the fluctuations at equilibrium (the thermal 
noise). Since it is clear that microreversibility is bro- 
ken beyond the linear response regime, it thus natural 
to ask whether higher-order fluctuation relations exist in 
the presence of an external field. We find an approxi- 
mate verification of such relations to next order in the 
voltage expansion. In other words, the asymmetric part 
of the second-order conductance and that of the linear- 
order noise are related to each order via a nontrivial fluc- 
tuation relation. Remarkably, our results are interesting 
because we treat interactions beyond the Hartree-Fock 
approximation. 

The paper is organized as follows. In Sec. [TT] we in- 
troduce the system and its Hamiltonian. Section. IIIII is 
devoted to the transport properties. We also analyze two 
approximations: the noninteracting case and the Hartree 
approach. Both limits have serious drawbacks and can 
even lead to unphysical predictions but we nevertheless 
include discuss them for pedagogical reasons. A detailed 
study of the transport coefficients in the Coulomb blo- 
cakde regime is contained in Sec. IIVI We examine the 
even-odd properties of the conductance terms both for 
zero and nonzero temperatures. In Sec. [V] we consider 
the Kondo regime, which is the relevant scenario at very 
low temperatures. We briefly discuss the limits of tem- 
peratures lower and higher than the Kondo temperature 
and then present numerical results for the magnetoasym- 
metric differential conductance as a function of the back- 
ground transmission and bias voltage. In Sec. IVII we 
calculate the noise and show that it is an asymmetric 
function of the magnetic flux. Finally, our results are 
summarized in Sec. IVIII 



Here, 



ct=L / R,k,cr 



E (We^c Rkl<7 c Lkr7 +h.c) , (4) 



k,k f ,a 



describes the two leads and the direct channel that cou- 
ples them. The gauge is chosen in such a way that an 
electron wave picks up the phase ip whenever it passes 
along the upper arm. The dot electrons obey, 



Kd = 22 e d d\d a + U n d1 -n dir . 



(5) 



while the tunneling Hamiltonian between the reservoirs 
and the dot is given by, 



h t = E y^ja + h, 

a=L/ R,k,a 



(6) 



where the tunneling amplitudes V a are asummed to be 
independent of k for simplicity. The same assumption is 
made for the direct transmission W. 



III. TRANSPORT PROPERTIES 

In the stationary limit, the current I can be calculated 
from the time evolution of the occupation number of the 
right contact (ur), 



I = h. 



dt 



[H,n R ]), 



(7) 



where n R = J2ka c Rka c Rkv- Using the Keldysh 
formalism, 49 the current becomes 



II. THEORETICAL MODEL 

The mesoscopic interferometer consists of an 
Aharonov-Bohm ring with a quantum dot inserted 
in one of its arms. The interferometer is coupled to 
left (L) and right (R) leads with continuous energy 
spectrum e k(T where k is the wavevector. An electron 
can travel either through the nonresonant arm with 
probability amplitude W or via the quantum dot with 
hopping terms V a where a — L, R is the contact index. 
The spin-degenerate level dot is denoted with e d and 
the charging energy is U. Finally, the magnetic flux $ 
piercing the ring results in an Aharonov-Bohm phase 
if = eQ/h. Hence, the Hamiltonian reads, 



H = H c +H d +Ht- 



(3) 



'■IE 



dui 



( V RG d a,RqA £ ) - V RGRqa,da(u) 



+ (We^g<^ Rq Jco) - W*e-^g<^ Lp Jco))} , (8) 

with the following definitions for the lesser Green fucn- 
tion(G<): 



(9a) 
(9b) 
(9c) 
(9d) 



In the case of energy-independent couplings or for pro- 
portionate couplings, the expresssion for the current is 
more conveniently recast in terms of a generalized trans- 
mission function Ta-(uj) for an electron with spin a and 
energy lu, 



G Lpa,Rqa 


= i 


\ C Rq<7 CL P<7 


G Rqa,Lpcr 


= i 


\ C Lpa C Rg (J 


^ da ,Rqcr 


= i 


\ C Rqad<r) > 


^ Rq<7,d<7 


= i 


[dtc Rqa ) . 



-Y 



duj %(u))[f L (ui) - /r(uj)} 



(10) 
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where and fn are the Fermi-Dirac distribution func- 
tions in the leads L and R, respectively. The transmission 
reads 



%(u) = T b + y/aT b n b cos(ip)TBe [S^M] 

-i{ a [i-r b cos 2 M] -r b }rim[a5 Mff M] . (n) 

Here, 71 is the transmission probability between the two 
leads along the direct channel, 



(i + 2 



(12) 



where £ = ir 2 W 2 plPr with Pl{r) the density of states for 
lead L(R). The reflection probability lZ b is determined 
from 7l b = 1 — %. The broadening of the dot level due 
to hybridization with states of lead L(R) reads Turn — 
2t|Vl(.r)| 2 /Ol(.r), where Pl(r) is the density of states for 
lead L{R). The total linewidth is V = T L + T R . In the 
wide band limit, we take Pur) = pa (and, consequently, 
£ and T) to be energy independent. In the presence of 
the upper bridge, the broadening becomes renormalized, 



r = 



r 



(13) 



Finally, the factor a = ATlTr/T 2 in Eq. (fTTj) quantifies 
the tunneling asymmetry (0 < a < 1). It yields a = 1 
for symmetric couplings (r^ = Yr = T/2). 

The problem has thus been reduced to the calculation 
of the dot retarded Green function Q r da da (u>)- For later 
convenience, we introduce the following definitions: 



({A(t),B(t')))< = 
((A(t),B(t')))> = 
((A(t),B(t'W = 
((A(t),B(t'W = 



-i{B(t')A(t)) , (14a) 

-i{A{t)B(t')) , (14b) 

-ie(t-t')([A(t),B(t')} + ), (14c) 

-iQ(t'-t)([A(t),B(t')]+), (14d) 



where r (a) stands for "retarded" ("advanced"). Using 
the operators A = d„ and B = d\ we obtain the dot 
Green function, 



dte l ^{{d(t),d\Q)))\ (15) 



where we have set t' — since the Hamiltonian is time 
independent. The equation of motion for Q r da dr7 is found 
as, 



1 i — — i ~ 

w - Ed + 7 V ®T b T cos(ip) + -T 



((d CT ,4))'- = i 



U((d M ,dt)) r . (16) 



puted, 



ed + iVS7lrcos(^) + |r 



(17) 



Substituting in Eq. (jllj) . the transmission probability be- 
comes, 



ti 2 + \v 2 



where 



u> = uj — Ed + -\/aT b T cos(ip) . 




'TZb cos(ip) + i sm(ip) 



(18) 

(19a) 
(19b) 



Equation (fl~8)) is evidently of the Fano type^ The Fano 
antiresonances arise as a consequence of interference be- 
tween a direct path channel (the upper arm in Fig. [T]) 
and a hopping path via a quasi-localized state (the quan- 
tum dot in the lower arm). As a result, a characteristic 
asymmetric transmission lineshape is obtained, which is 
described with the (generally complex) Fano parameter 
qf- When ip is a multiple of it, the transmission vanishes 
at the special energy point given by e(e) = —qpT/2. On 
the other hand, for vanishingly small transmission along 
the direct channel (Tb — > 0), Equation (fT5)) reduces to the 
Lorentzian form of the transmission resonance through a 
noninteracting dot. 

We take a bias V symmetrically applied to the elec- 
trodes (fj,L — ~pr = eV/2) and insert Eq. (fT5)l in 
Eq. (|10l) . Next, we expand the current -voltage character- 



istics, / = GqV + GiV 2 + G 2 V 3 + ■■■ = £„> G n V n+1 . 
We find that the even coefficients Cm (n — 0,1,2...) 
are functions of cos ip while the odd coefficients vanish, 
G2n+i =0. As a consequence, the current is an even 
function of the flux and fulfills the Onsager symmetry. 
This result also holds for finite temperatures. For in- 
stance, the linear conductance reads, 

2e 2 /3 



Go = [v^aTbTlb cos(^)Im [* ] 

11 Z7T I. 



- [a (1 - % cos 2 fo>)) - T b ] Re [tf „] \ , (20) 
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where j3 is the inverse temperature and 
* (l + + ^Sr) denotes the digamma function. 
Note that we have subtracted the offset term 
f dui 7&[/i,(w) — fn(uj)]. In the zero temperature case, 
Go becomes, 



A. Noninteracting case 

In the absence of the interaction, one sets U = in 
Eq. (fT6"|) and the retarded Green function is readily com- 




TsdVoiTbllb cos(y>) 



£2 

4 



(Sd) 2 

f 2 [a (1 -% cos 2 -%] 



(Sd) 2 + £ 



(21) 
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where 



£d — £<J 



'aTbT cos(y) 



(22) 



Equation (|2~T|) is clearly an even function of (p. These 
results show that in the absence of interactions, trans- 
port is B-symmetric to all orders in voltage. However, 
a word of caution is in order. Neglecting interactions in 
the nonlinear regime of transport can lead to unphys- 
ical results (e.g., gauge invariance can be broken) i 47 i 50 
Therefore, to give reliable results away from equilibrium 
we must include interactions at least in the lowest level 
of approximation. 



B. Hartree approximation and even-odd behavior 

We now introduce interactions in the most simple way, 
namely, we use in Eq. (|16|) the following decoupling, 



{{d a n ds ,dl)) w (n d9 ){{d a ,d a )) 



(23) 



This Hartree approximation is well known to sponta- 
neously generate local moment formation^ in the quan- 
tum dot. Although this result is physically meaningless 
since the Hamiltonian [Eq. is invariant under spin ro- 
tations, the approximation is useful as a benchmark for 
more elaborate models. 

The retarded Green function is found to be, 



'da.da 



1 



u-e d - U{n ds ) + IVTbT cos(^) + §r ' 



which was anticipated in the Introduction. The electron 
occupation (n da ) can be obtained from 



{ n do) = ^~ J duj Q daAa {uj) 



(25) 



In general, for interacting systems, the lesser Green's 
function cannot be directly obtained from the equation- 
of-motion technique without introducing additional as- 
sumptions. However, we note that only the integral of 
G dada {u) is, in fact, needed in Eq. (|25|) . This observa- 
tion allows us to bypass any approximation involved in 
computing G d<rM (u), yielding^ 



\n da ) = 



1 



(26) 

An alternative derivation is presented in App. A. In 
Eq. (|26[) . f peq (io) denotes a pseudoequilibrium distribu- 
tion function which is not, quite generally, of the Fermi- 
Dirac type, 



Setting Q a da4 M = we see that E( l s - <©> 

(|26| and (|27| form a closed system of equations which 
must be solved self-consistently. But before proceeding 
with such a calculation, we point out to the presence of 
a tp- asymmetric term already in Eq. (|2T[) . Note that this 
term is nonzero only in the nonequilibrium case (/l ^ 

/*)■ 



1. Zero temperature case 

For = T r, the pseudoequilibrium distribution func- 
tion can be simplified, 



[f L (u) + /rH] + V7Tsin(y) [f L (u) - f R (u)} 

2 

(28) 

and from Eq. (|26[) we find an exact expression for the 
occupation, 



(n do 



— {[l + y/%Bm(<pj\) 

2 (e d + U{n ds ) - /j, L ) 



x cot 



1 - V 7h sin(<p) 



cot 



2 (e d + U (n ds ) - fi R ) 



(29) 



As introduced in Sec. HI injection from the left lead con- 
tributes to the dot occupation with a term 1 + y/Tb sin(ip) 
while the contribution from the right lead is given by 
1 — v7f> sin(^). Both terms cancel out at equilibrum, re- 
gardless of interaction, but survive in the presence of a 
finite bias. 

We solve Eq. (f2"9"]l iteratively. We write, 



(ndo 



(ndc 



\(o) 



(ndc 



(eV) 



(nda 



,(2) 



(eVf 



(30) 

insert this expansion in Eq. (|29p and assume (n da ) = 
(n d g). Then, we obtain the following expansion coeffi- 
cients, 



(rw) (0) = 



1 

— tan 

7T 



2? ! 



{ndo) 



(i) _ 



2{- [(^f + ?]+|^} 



(31a) 
(31b) 



(ndc 



v(2) _ 



r^[i + 4(f/(^)W) 


2 - 4Uy/%sin(ip)(n a 




4 |tt 


\nf + ?; 




(3- 


Lc) 



f P eq(uj) = h j^fy / a^sin((/?) [f L (uj) - f R (u)] 

- j^-^ [{r L + tr R ) f L ( u ) + (r R + ^r L ) /«( w )] 



(27) 
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(n d<7 )& =3!|2f[/ 
+3(n d<r )W [l2(^) 2 - 



4U' ( 12(^) 2 - r 

f 2 



((^> (1) ) : 



2 



r(-12(^) 2 +r 2 ) (l + 12(tf<n*,)«) 
f 2 ) (n dCT )( 2 )]} 



-12rf/g° (4(eS) 




x< 192 



where we have defined 



(31d) 



r 6 r C os(^), (32) 



which is an even function of the field. ^From Eq. pip , 
we infer the relations, 



(rw) (2) M = +(rw) (2) (-^) 



(33a) 
(33b) 
(33c) 
(33d) 



Therefore, we can then conclude that the nonequilibrium 
charge response of the system is not a symmetric function 
of the flux, as shown in the odd coefficients in Eq. (|33bj) 
and Eq. (|33d|) . In fact, these coefficients are antisym- 
metric when the field is inverted. On the other hand, 
the even coefficients are all symmetric under reversal of 
ip. This restoration of the Onsager symmetry for the 
even coefficients of a current-voltage expansion has been 
observed experimentally in rings (see Ref. I25I ). In this 
section, we do not attempt to make a direct compari- 
son with the experimental data since the conductance 
within the Hartree approximation would give wrong re- 
sults due to the aforementioned breaking of the spin ro- 
tation symmetry. Nevertheless, in the next section we 
demonstrate that this even-odd behavior is not an arte- 
fact of the Hartree approximation and persists in a better 
treatment of Coulomb interaction from which a physi- 
cally meaningful conductance can be extracted. 



2. Nonzero temperature case 

For completeness, we briefly discuss the nonzero tem- 
perature case. Equation (|29|) is replaced with 






1 






X 




— Im 


H 




2 ~ 


7T 





We now substitute the expansion Eq. ([30|) in Eq. 
and find exactly the same relations as Eq. (I33b[) . 
instance, the first two expansion coefficients read, 



\rida 



(Uda 



,(0) 



1 

2^ 



1 



Im [* ] 



a y/%sm(tp)Re 




4tt3 1 + |^ Re 







(34) 
For 



(35) 
(36) 



(n) 

where ^ = 



(l + f£ + ^H) is the polygamma 
function (the nth derivative of the digamma function de- 
fined above) 4£ We again see the antisymmetric charge re- 
sponse of the system due to the sin ip term in the leading- 
order nonequilibrium coefficient (rida)^ ■ 

IV. COULOMB BLOCKADE REGIME 

In the Coulomb blockade regime of two-terminal quan- 
tum dots, transport takes place only through two reso- 
nances approximately located at e d and e d + U . Clearly, 
the retarded Green function given by Eq. (j2"4")) does not 
show this behavior and consequently we must perform a 
higher-order truncation in Eq. (|16[) . This way, one ob- 
tains the equation of motion for ((d a nda,d\)): 

{ui-e d - U) ((d a n d a, d\)) = (n ds ) 



a.k 



a , k 



E^ c lfc ff ^4))- (37) 



a.k 



To obtain the two-peak solution, we keep only the 
first term on the right-hand side of Eq. (j3"T)l , calcu- 
late its equation of motion, and make the following 
approximations 



((Cakacl, k , s d s ,dl)} « 0, 
{{CctkcrdgCa'k'a,^)) ~ 0. 



(38a) 
(38b) 
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Then, 



A further expansion of / in powers of V finally gives, 



1 - (rtda) 



{nds) 



G 2 n(^) = G 2n (-f) , 
G2n+l{f) = -G 2n +l(-(fi) ■ 



(47a) 
(47b) 



e d -U + lV^nrcos(ip) + iT 



(39) 



We note that the retarded Green function now correctly 
shows two peaks located at e d — jy/aJbT cos(ip) and e d — 
-jV 'a7&r cos(</j) + f7 with weights 1 — {n ds ) and (rids), 
respectively. 

Using Eq. (|26p we find that the occupation is given by 



{rida) = 



4 + r L {X L (go) _ )) + Tfl (x a (gO) _ )) 

(40) 

where we have used the following definitions, 



T L {K) 

X a (x) 



(i + Or 



± y/ a% s'm(ip) , (41a) 



1 Im 

7T 



2 + Air + 2tt 



t rl 



£d - | aTbT cos((p) , 

£d + 17 - -7 y/a%T cos(ip) , 



(41b) 
(41c) 

(41d) 



For symmetric couplings (r^ = Yr) or a completely 
open nonresonant channel (£ = 1), we find for the par- 
ticular case of symmetric bias (/il = — fi R = V/2) that 



(n da )(ip,V) = (n<ur)(-<P,-V) 



(42) 



Physically, this corresponds to an invariance of the whole 
system when both the magnetic field and the electric bias 
are inverted^ A series expansion in V then yields, 



(nd„) {0 Hv) + M {1) (v)V + --- = 

(n dlT )^{-ip)-{n dcr )^{-<p)V + - 

Thus, we have 



(43) 

(44a) 
(44b) 



These equations represent a generalization of the Hartree 
case [Eqs. (13"3"|) ] to the Coulomb blockade regime. It then 
follows that 



and 



%(<p,V) = %(-<p,-V) 
I(<p,V) = -I(-<p,-V). 



(45) 
(46) 



i.e., the even (odd) conductance coefficients are symmet- 
ric (antisymmetric) functions of the flux. We see here a 
crucial difference compared to the noninteracting case 
discussed earlier. For U = the current is always a sym- 
metric function of tp regardless of the applied voltage. In 
the interacting case, we find that the odd coefficients of 
the conductance are not invariant when the field orienta- 
tion is inverted. 

The above property is not general and can be traced 
back to the spatial symmetry of the system (symmetric 
couplings and symmmetric bias). In the more general 
case of asymmetric couplings (r^ ^ T^), the occupa- 
tion symmetry given by Eq. (j42|) is not fulfilled. To see 
this, for simplicity we take the limit U — >■ oo. Then, the 
occupation is given by 



{n d c) 



X L +X R + A(<p) {Xl-Ir) 



A + X L +X R + [s/n b ST + ^/W b sin(p\{I L -T R ) ' 
where X a is evaluated at x — eS and we have defined 



A(ip) = y/TZ b 5T + yf<W b sm(<p) 



and 



ST = 



Tl-Ti 

r 



(49) 



(50) 



We now expand the occupation as a function of V and 
find that the expansion coefficients are given by 



{ n a 2n) = Y,ct\ V )A^) 



(51a) 



»./n--- 1> = E C !Ki 1 W am+1 (v), (51b) 



m=0 
m=0 

with the C's fulfilling, 



(52) 



We give the explicit expressions for the first leading-order 
coefficients, 



7r - 2 Im * 



ii 



{rider 



{ride 



\( 2 ) 



3?r - 2Im* 

(37r-2Im* ) 2 //3 ' 

(3?r - 2Im* )Im*[ ) 2) -4 



Re* 



(i) 



(53a) 
(53b) 



87r(37r -2Im* ) IP 2 



(53c) 



Now, the dimensionless function A(ip) is a small quan- 
tity for almost all cases. We show in Fig. [2] that A is 
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FIG. 2. (Color online) A(<p) as a function of 71 and ST for 
ip = Tv/2. Refer to Eq. |g9j) . 



always smaller than 1 even for the case <p = tt/2. As a 
result, we can safely neglect A™(<^) for all n > 1, thus 
keeping the first-order term only. This implies that the 
even expansion coefficients are always symmetric under 
the reversal of if but the odd ones do not show any par- 
ticular symmetry since Tl ^^r- Since we know that the 
transmission 7^-(w) from Eq. (|10[) obeys the same sym- 
metry as the occupation (ri^o-), it follows that 



G 2n (ip) = G2n(-<P) , 

G 2n +i(f) ^ -G 2n +i(-<p) ■ 



(54a) 
(54b) 



This is precisely the behavior that was observed in Ref.l25l 
for an asymmetric ring. In Fig. [3] we numerically cal- 
culate the first four conductance coefficients in the cur- 
rent expansion for two different values of the background 
transmission. In both cases, Go obeys reciprocity, as ex- 
pected. For a partially open direct channel (Tb = 0.5), 
the leading-order nonlinearity, G\, is magnetoasymmet- 
ric but the Onsager symmetry is recovered for G 2 and 
later destroyed again in G3. These results are in agree- 
ment with the experiment^ In the case of a fully open 
direct channel (Tb — 1) the odd coefficients are still asym- 
metric but they are now odd functions of the magnetic 
flux, in agreement with Eq. (|47b[) . 



V. KONDO CORRRELATIONS 

We can now go to next order in the equation-of-motion 
technique to describe the onset of Kondo correlations. 
Thus, we obtain the equations of motion for the three 
functions appearing on the right-hand side of Eq. (|37|) . 
and approximate the new Green's functions that appear 
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FIG. 3. (Color online) Conductance coefficients as a function 
of ip. Go and G2 are offset for simplicity. Parameters are To = 
1, ST = 0.2, r Q = r (l ± ST), e d = -5r , and T = r . Here, 
the conductance coefficients are scaled by e" +2 / ' (2 n ~ 1 n\h). 



in the procedure by making the decouplings 
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({c ak ,,c\ rii dv,dl)) w (cl r9 d s )((c akt7 ,dl)} , 

({Cakadycp r a,dl.)) PS (d\c [jrS ) ((c akcr , d\)) , 
{{d\ c akaCpr<i,dX)) W (d\c ak d) ((cpra, , 

(4 rff c Qfcff )((d <7 ,4)) , 

)«d*,4)), 
d 9 ){{cp r<T ,dl)) ■ 



({c^ rff Cak 9 d a ,dl)) 

(( C lk* C Pr*d ai dD) 
{{ c ak<xd»C0r<T,dl)) 



(55a) 
(55b) 
(55c) 
(55d) 
(55e) 
(55f) 



In what follows, we take the limit U — > 00, in which 
case the term ((c i aks d r7 da,d'l)) does not give any contri- 
bution. After little algebra, we find, 



1 - (rids) - 5n ds (uj) 



where 



uj + i0+ — ed — (1 — 8n ds (ui)) Eo(w) - 

(56) 



5n ds (u) = -— dui 7 

Z7T / U>' 



(57) 
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and 



S (w) 
Si(w) 



r 1 r— . . 

~ l 2 ~ aibF cos(tp) , 

f 

"2^ 



(58a) 



, , fpeq(u') 

duj w /_ u ,_ i0 + 



{l+ [Eo(a/)«£fe,4»^]*} • (58b) 



The derivation of this expression for the retarded Green 
function is explained in App. |B] 

To lowest order in T, Eq. (|56[) can be further simplified 

as 



1 - (n-da) 



10 + i0 + — Ed — T, r (uj) 



where 



with 



E r (io) = So(w) + Si(w) , 



fpeqjw') 

-w-iO+ 



(59) 



(60) 



(61) 



Here, we note that for = the self-energy obeys the 
following symmetry 



The Kondo temperature marks the energy scale below 
which nontrivial spin fluctuations start to play a domi- 
nant role, leading to an antiferromagnetic exchange be- 
tween the dot electron and the conduction electrons. We 
note that for a quantum dot inserted in an Aharonov- 
Bohm ring and in the presence of an applied flux, Tk 
depends on £ and ip but the dependence on the flux is 
weak^ Furthermore, the Kondo temperature is a static 
quantity and, as such, is always a symmetric function of 
the flux. 

We compare in Fig. 2] the value of Tk with the Kondo 
temperature of a two-terminal quantum dot, 



k R T, 



(0) 



BJ- K 



n / 2lT£d X 
exp I ~Y~ J 



(66) 



and plot the ratio Tk/T^ as a function of the direct 
channel tunneling probability [Fig. Ufa)] and the flux 
[Fig. HJb)]. For tp — and % = we recover Tk = 
T^ as expected. As Tb increases for ip = 0, the Kondo 
temperature decreases since electrons preferably travel 
along the upper arm. However, for fluxes above tp — 
7r/2, the level renormalization due to cos 93 is positive 

and the curve Tk/T^ becomes nonmonotonous due to 
the competition between the renormalized dot level and 
broadening. This can be more clearly seen in Fig. HJb). 



E r (^V) = £ r (-^-V). 



(62) 



In turn, this property implies that the occupation and 
the conductance obey the even-odd symmetry also in the 
Kondo regime (at least when the coupling is not very 
strong). 

Together with Eqs. (|25l) and (|27p the Green function 
of Eq. ([56]) can be obtained from a self-consistent proce- 
dure. But before solving this system of equations using 
numerical methods, we briefly discuss two limits (high 
and low temperatures) to clarify the origin of magne- 
toasymmetries in the Kondo regime. 



A. High-temperature regime 



B. Low-temperature regime 

At low-temperature 5n ds must be large, especially near 
the Fermi level. If we suppose that Q r da dcr (w) varies 
smoothly near the Fermi level, Eq. (|57|) can be approxi- 
mated as 



5n ds {ui) rj -— [G r d a,d„]' 



in , 2tt 
h In 

2 (3D 



Inserting Eq. (|6"7| into Eq. (|56p . we find 



(67) 



In this case, 8nds is a small correction and an expan- 
sion can be done. To first order in T it can be shown that 
the position of the virtual level is renormalized to e' d 



e d = e d --s/aTbTcos{ip)-—\n^ , 

(63) 

for (}(u) - E F ) > 1. If (i\e d - E F \ > 1 we have 

E F \ 



I y \e 

e' d = e d - -VaTbT cosftp) - — In 



(64) 



From the equation above, the Kondo temperature is 
given by 



knT, 



D exp 



2?r / 1 ,— . . 
— \ £d ~ 7V a lbT cos(yj) 

r V 4 



(65) 



(68) 



Here, the value of 9 is related to the number of d electrons 
according to the Friedel-Langreth sum rule^ 



TT{nd) 



7T 

2 



This implies 



((d a ,dt)) Ep = J-. 

II I A 

Then, the linear conductance can be written as 
G( >M = ^a[l-T h cos 2 M] . 



(69) 



(70) 



(71) 
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FIG. 4. (Color online) Kondo temperature Tk as a function 
of 71 and tp. The parameters are D = 1, Ed = —0.05, U = oo, 
and T L =r R = 0.031. For these parameters, T { ° ] w 0.0063. 
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FIG. 5. (Color online) Differential conductance dl/dV 
(2e 2 /h) as a function of the applied bias eV for D = 1, 
e d = -0.05, U = oo, T L = F R = 0.031, and tp = tt/4. 
Temperatures are T = 10T^ 0) (a) and 0.1T) f 0) (b). 



This result is exact in the limit ksT^V — > 0. For nonzero 
voltages, one should take into account that the imaginary 
part of the interaction self-energy of Q r drj dcr depends on V 
but this dependence is weak for eV <C IzbTk and can be 
safely neglected. As a result, deep in the Kondo regime 
the conductance preserves the Onsager symmetry since in 
the Fermi liquid picture the Kondo resonance behaves as 
a noninteracting system with renormalized parameters. 
Charge fluctuations are quenched and transport becomes 
B-symmetric. This regime is beyond the scope of our 
method and we prefer not to present numerical results 
for very low temperatures. However, the expected sce- 
nario would be as follows: for very low temperatures the 
current would be S-symmetric and asymmetries would 
arise as temperature approaches T% . In the opposite 

case, for temperatures much larger than tJ^ transport 
is thermally assisted and the magnetoasymmetric effect 
also disappears^ Therefore, we expect a large magne- 
toasymmetry for temperatures of the order of for 
which charge fluctuations are large. We confirm this ex- 



pectation in the numerical results reported below. 

C. Numerical results 

We numerically investigate the evolution of the mag- 
netoasymmetry when temperature is lowered from the 
Coulomb blockade regime to the Kondo temperature. We 
first illustrate the generic behavior in Fig.[5j which shows 
the differential condutance dl/dV a a function of the ap- 
plied voltage V for opposite orientations of the magnetic 
field. To compute the derivative of the current we have 
employed a numerical finite difference method. 

In the top panel of Fig. [5j the temperature is large 
enough that Kondo correlations can be neglected. Then, 
the dot is in the Coulomb blockade regime and a small 
current is expected since the dot level is below the Fermi 
energy (sd = —0.05). However, the bridge channel is par- 
tially open and the system conductance reaches around 
0.55 x 2e 2 /h at V = 0. This value is independent of the 
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magnetic orientation, as expected. But when V departs 
from equilibrium, the differential conductance behaves 
differently for +ip and —(p. We recall that the effective 
position of the effective resonance depends on the charge 
state of the dot, as discussed in Sec. IIV1 Since the charge 
response of the system is not a symmetric function of (p, 
dl/dV peaks at different voltages for opposite magnetic 
fields. 

In Fig. EJb) we depict dl/dV in the low temperature 
case. We observe for both field orientations a dip around 
V = 0. This dip is known to arise from the destructive in- 
terference between partial waves propagating through the 
upper arm and resonantly hopping across the dot^ We 
emphasize that the dot bare level, Ed, is the same for both 
calculations but in the Kondo regime transport is domi- 
nated by the narrow resonance formed at the Fermi level 
due to the higher-order tunneling processes that origi- 
nate the Kondo effect. As a result, the Fano interference 
between the Kondo resonance and the background chan- 
nel gives rise to the pronounced dip at zero bias. In our 
case, we obtain an asymmetric lineshape for the dip due 
to the magnetoasymmetric response of the dot away from 
equilibrium. As a consequence, the difference in dl/dV 
for +ip and — <p is more visible in the Kondo regime, as 
can be seen in Fig. E2h) compared to Fig. EJa). 

We now define the symmetric (+) and antisymmetric 
parts (— ) of dl/dV, 



0.04 



dl± _ dljip) dl(-ip) 
W ~ dV dV 



(72) 



In Fig.|6]we plot the ratio between these two components 
as a function of the applied bias V for a fixed value of 
the flux (tp — 7r/4) and for different values of the nonres- 
onant transmission %■ In Fig^a) we set the tempera- 
ture to a high value compared to the Kondo temperature. 
For Tb = 0.5 the magnetoasymmetry is always finite for 
V 0. For voltages around zero, the magnetoasymme- 
try is a linear function of V since the largest contribution 
stems from the G\ coefficient in the current -voltage ex- 
pansion. In the limit of high bias, the magnetoasymme- 
try saturates?^ Interestingly, with increasing Tb the mag- 
netoasymmetry is reduced and changes sign for a fixed 
V. Thefore, the sign of the asymmetry can be tuned with 
the background transmission of the nonresonant channel. 
The situation is similar to the magnetoasymmetry of a 
two-terminal quantum dot when transport is dominated 
by elastic cotunneling processes^ In that case, the sign 
of the asymmetry can be changed with the gate voltage 
which moves the dot level position above and below the 
particle-hole symmetric points In our case, Tb acts as 
an effective gate which changes the position of the level 
since Sq is renormalized according to Eq. (|22[) . 

When temperature is lowered, we observe that the 
transition from positive to negative asymmetries as V 
is tuned, is rather abrupt, see Figjfjb). We note that as 
voltage approaches V — one sweeps along the strongly 
asymmetric dip structure found in Fig E3h), for which 
the difference between the cases +(p and —tp is most clear. 
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FIG. 6. (Color online) Ratio of differential conductances 
dI±/dV (2e 2 /h) as a function of the applied bias eV . Here, 
dI±/dV is denned as d(I(<p) ± I( — ip))/dV. Parameters are 
D = 1, e d = -0.05, U = oo, F L = r fl = 0.031, and (p = ir/4. 



Temperatures are T = 10T^ 0) (a) and 0.1T£> (b) 



n(0) 



Then, it is the Kondo resonance that produces the abrupt 
change in the magnetoasymmetry profile as compared to 
the high temperature case [FigG^a)]. As a consequence, 
Kondo correlations enhance the deviations from the On- 
sager symmetry since the narrow resonance is more sensi- 
tive to changes in the orientation of the magnetic fields 
For instance, we observe in FigG^b) a revival of the mag- 
netoasymmetry for Tb = 0.75, which almost vanished in 
the high temperature case. However, if temperature is 
further lowered (T -C Tjg) for Ed -C Ep, charge fluctua- 
tions would be quenched and the Kondo resonance would 
be pinned at the Fermi level, independently of <p and Tb- 
As a consequence, the magnetoasymmetry would tend to 
vanish. 
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VI. SHOT NOISE 

The shot noise is a valuable tool in the characteriza- 
tion of the transport properties of mesoscopic systems^ 
For systems described with Anderson impurity models 
like ours, the electron repulsion term U introduces cor- 
relations which can be investigated through the noise. 
Then, the problem becomes involved, although the effect 
of Kondo correlations in the shot noise have been already 
addressed in a number of papers^ - — . 

Magnetoasymmetries in noise have recently attracted 
a good deal of attention due to the (weakly) nonequi- 
librium relations between the asymmetries of the cur- 
rent and that of the noise to leading order in a voltage 
expansion!^ - — The subject is also of interest because 
it poses questions about the validity of fluctuation theo- 
rems out of equilibrium^ In this section our goal is to 
calculate the noise power for our system in the limits of 
both weak and strong electron-electron interactions and 
check the nonequilibrium fluctuation relations. 

The current noise between terminals a and j3 is defined 

as 

s a ^t-t') = ^{{[i a (t)J (t')} + )-2{i a }{i p }] , (73) 

where I represents a current operator. The Fourier trans- 
formation of the current noise (the noise power) reads 



S(e>) 



dt 



(74) 



where we have defined the Fourier transform without the 
prefactor 2. In the following, we present results for the 
zero frequency case [the shot noise 5 = 5(0)]. 

The noise definition of Eq. ([73"]) contains correlations 
between currents that, quite generally, involve four oper- 
ators. To treat the resulting two-body Green's functions, 
we make use of cluster expansion ; 68 ' 70 ! 74 
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linear regime of the noise response. In Fig. [7] we show 
the noise as a function of V in the strongly interacting 
case. We can observe that the slope of the noise curves 
at V = differ for opposite field orientations. Another 
interesting feature is that for some voltages the nonequi- 
librium noise can be reduced from its equilibrium valued 

To gain further insight, we expand the noise in powers 
of V, 



S = Sq + S!V + 



(77) 



5o is the equilibrium noise describing thermal fluctua- 
tions. Since these fluctuations do not distinguish between 
+ip and —ip, Sq is an even function of the magnetic field. 
An alternate proof of this statement is based on the equi- 
librium fluctuation-dissipation theorem, which relates Sq 
to the linear conductance Go, 



+ <W<o 



,>*')(d 



o„ 



fi'cr' I 



(75) 



which amounts to neglecting two-body connected Green's 
functions. As a result, the shot noise is expressed in 
terms of one-body Green's functions. This is a strong 
assumption that can lead to deviations from well known 
relations. Nevertheless, interactions are included (at the 
level of the Lacroix's approximation). The calculation 
is lengthy and we refer the reader to App. C. Here, we 
consider limit cases only. 

For the noninteracting case and at zero temperature 
we recover the known expression, 



"E [ eV2 ds%(u) [l -%(")} , (76) 

1 a J-zVjl 



where the transmission is given by Eq. (j!8[) . As expected, 
the noise is an even function of ip to all orders in V. How- 
ever, interactions destroy this symmetry already in the 



So = 2k B TG • 



(78) 



Since for Go the Onsager symmetry holds, So should be 
even for a two-terminal setup. 

We now show the numerical results of our model for 
different temperatures. We define the symmetric and an- 
tisymmetric components of the noise and the current as 
before, 



5± = S(+<p) ± S(-p) , 
7± = I(+p) ± I(-<p) . 



(79a) 
(79b) 



We consider symmetric couplings. As result, the even- 
odd properties of the transport coefficients allow us to 
write, 



5. 



I+/V 

5_ 



= 2k B T 
= 2k B T. 



(80a) 
(80b) 
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dissipation theorem as a function of if. Parameters are D = 1, 
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Corrections to these relations are of order V 2 and can be 
neglected for eV <C feT. The first relation is merely a 
restatement of Eq. ([75)1 . The second relation is a nonequi- 
librium relation that connect the magnetoasymmetries 
corresponding to both the linear- response noise (Si) and 
the leading-order nonlinear conductance (Gi). 

We take a small voltage (V = 0.01T^ 0) /e) and plot in 
Fig. [5] the relation given by Eq. (|80al) as a function of the 
flux for two different values of the background transmis- 
sion. In bot cases we find that Eq. (I80a[) is approximately 
fulfilled with small deviations which we attribute to the 
assumption of Eq. (jT5j) . We note that deviations grow for 
smaller temperatures since in this case the Kondo corre- 
lations become more relevant and our model for the noise 
starts to break down. When the nonresonant channel is 
fully open (bottom panel of Fig. [5]), the deviations are 
less important since electrons preferably travel along the 
upper arm and consequently feel less the intradot inter- 
actions. 



The validity of the nonequilibrium fluctuation relation 
[Eq. (|80b[) ] is analyzed in Fig. [9j Here we plot separately 
the two terms of Eq. (|80b[) . We find a strong ressem- 
blance between 5i and G\ for all magnetic fields. Devi- 
ations also exist as in the calculation of the symmetric 
components but they fulfill the same pattern, namely, 
they tend to disappear when the background transmis- 
sion is close to 1 and the temperature increases well above 
the Kondo temperature. Although our results are not a 
conclusive proof of Eq. (180b)) for strongly interacting sys- 
tems, the errors are small and compatible with the same 
deviations found in the equilibrium case (Fig. [S]). 

VII. CONCLUSIONS 

We have shown that the current -voltage characteristics 
of a two-terminal quantum-dot mesoscopic interferome- 
ter is not an even function of the applied flux in the non- 
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linear regime of transport and when intradot interactions 
are taken into account. The interference pattern of the 
ring at a finite bias is thus not symmetric under reversal 
of the magnetic field. We have carefully investigated the 
symmetry properties of the conductance coefficients in a 
current-voltage expansion. Our discussions are based on 
the properties of the charge response of the dot when a 
finite bias is applied to the system. When the quantum 
dot is in the Coulomb-blockade regime, we find for most 
cases that the even coefficients are symmetric functions 
of the field while the odd coefficients does not show any 
relevant symmetry in the general case. Only when the 
dot is symmetrically coupled to the leads the odd coeffi- 
cients are antisymmetric. 

We have also calculated the magnetoasymmetry of the 
system in the strong coupling regime, when the dot is de- 
scribed with Kondo correlations. In this case, the mag- 
netoasymmetry shows an abrupt transition between pos- 
itive and negative values when the voltage crosses the 
Fermi energy. As a result, the Kondo resonance dom- 
inates the magnetoasymmetry lineshape when the volt- 
age is of the order of the Kondo temperature. A further 
extension of this work could be focused on the very low 
temperature regime using, e.g., slave-boson techniques. 

Finally, we have investigated the asymmetry in the 
shot noise, finding a correlation between the noise and 
the current magnetoasymmetry to leading order in the 
applied voltage. This nonequilibrium fluctuation relation 
seems to apply in a wide range of parameters (temper- 
ature, direct channel transmission and applied fluxes). 
However, further work is needed to reduce the deviations 
which arc most probably due to our approximations. In 
particular, it would be interesting to analyze the role of 
the third cumulant of the current (or better the entire 
full counting statistics) under field reversal. 
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Appendix A: Lesser Green's Function of the Dot 

In general, the lesser Green's function for an interact- 
ing dot cannot be obtained from the equation of mo- 
tion technique without introducing additional assump- 
tions. Following Ng's heuristic approach, here we em- 
ploy an ansatz for interacting lesser and greater Green's 
functions^. Using the lesser and greater self-energies for 



a noninteracting dot 

S<(°)H = +*^f Q / Q H, (Ala) 

a 

E>< >H = -iX;fa[l-/ o M], (Alb) 



with 



{1 l^ }2 (r L +gr fl ) + ^v^sin(^, (A2a) 

f «= (T^j2 (r^ + M-I-v/a^sin^). ( A2b ) 

we assume that lesser and greater Green's functions for 
an interacting dot can be written in the form 



£>(w) 



:^r a R*(u). 



(A3a) 
(A3b) 



The explicit form of i? cr (oj) can be obtained from the 
relation 

E>( W )-E< («) = £»-£•(«). (A4) 
Employing the identity, 

GLM - Gd*,dM = - Q%r,M > (A5) 

it finally yields 
with 



fpeq (w) 



(A7) 



Appendix B: Evaluation of expectation values 

In deriving Eq. (|56|) . we have to evaluate the expec- 
tation values {d\c aks ) and J2/3,r V ^( c lra c ak<j)- First, let 
us concentrate on {d\c aks ). In equilibirum, the quanti- 
ties can be calculated by using the fluctuation-dissipation 
theorem, 

{d\c aks ) = -- [ duj f eq (ui)lm ((c ak a,dl))l . (Bl) 



However, the fluctuation-dissipation theorem cannot be 
employed in nonequilibrium. Then, in general, 



{d\c aks ) = i fdu {{c aks ,d\))< . (B2) 
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Note that ((c aks , 4))w cannot be obtained exactly be- 
cuase of the dot lesser Green's function. Here, we thus 
assume the following pseudoequilibrium form 



E 



{{c akSl dl))< = -f peq (uj) («C Q ,fcg.,4»L - «Cafcg-,4» 



t\\a 



X 



uj + i0+ - £ Rks 1 + £ 



duo 



I fpeqiw') 



UJ' — UJ — I 



i0+ 



Then, we have 



(B3) 



2 ~ I duJ iWgM ((( C c*r>4}}L(</>) ~ ((c a fcS,4»°M 



where 



(vk + iVfe+^Vh) {l+ [Eo(^)((^,4)>L'] } • 

(B9b) 



Appendix C: Expression of the shot noise 

The current noise is defined as 



(B4) 



«cl**,4»: = ff| (Vl - iVte- zv V R ) {{d- a ,d\))l , 



(B5a) 



(B5b) 



and 



4»S = ffl (vk + «<fe,4»£ . 

(B6b) 

with g^g = i *0 + — £ ate)- Due to the approxi- 
mations, however, the retarded Green's function is not 
the complex conjugate of the advanced one so that this 
equation gives an unphysical result. To resolve this diffi- 
culty, we thus replace the quantity in the parenthesis by 
the time-reversal pair. That is, 

{d\c akS ) = J duJ fpeq{u) 

x (((c ak9 ,4))l(-uj) - ({c ak *,4)Z(v)) • (B7) 



This yields a real expectation value as should be. Using 
Eq. (|B7|) . we thus have 



E 



(dic Lk a) 



Po 



duo 



I /peg(w') 



UJ + i0+ - ELkS l+U Uj'-UJ-iO+ 



E 



{d\c Rks ) 



= Po 

uj + i0+ - e Rks 1 + £ 



duj 



I fpeqjuj') 

uj' -uj- i0+ 



x (y R + i^le +i *V L ) \((d-„,dl))l] ■ (B8b) 
In the same way: 



E 



PO 



uj + iO+- e Lks 1 + £ 
. rz _„•,„„ \ r. r„ 



duj 



uj' -uj- i0+ 

, A , , 1*1 



S^(t-t') = ^{<[4(t),^(0]+)-2<J Q )(^)} . (CI) 



The current operator can be calculated from the time 
evolution of the occupation number operator 



f E [Vsc^da - V^4c Rqa ] I . (C2) 

qG-R,cr I 



The current can be then obtained by taking the average 
(B6a) of I R . Using the Keldysh G reen's functions, the current 
can be expressed as 



-We^G< q ^ Lpa (t,t)) 

+ E (*,*)" VrO^M (*>*)) • ( C3 ) 

q€R,<J 

The Fourier transformation of the current noise is 



(C4) 



The zero-frequency noise power S = 5(0) is referred to 
as shot noise. In the derivation, we make use of cluster 
expansion to treat two-body Green's functionsi 68 i 70 ' 74 

+ <W<6t 6 vl<JI ){d v<T 6\ l(; ,) (C5) 



In the frequency domain, using Eqs. (|C2p and (|C3[) and 
cluster expansion the shot noise is given by 
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S = 



1 2(ief 



2tt 



duj {W 



2 '^Gt^M^ 



Lp(T,Rq' a 



1 Rq' a,Lp<r 



( w )£fl 9 <x,Lp' CT M] 



W2 [Gi q >a,Rq<y( UJ ) Q Lpa,Lp>A UJ )+Qtp><j,L P <y 



W 
W 
W 
W 



V^g< q , atRqa (u;)GL P a4^) + yR^R q :,R q MQlr,L P 'M 

v^g< ql ^ Lpa {u)g> qaAa {uj) + v R e + ^g< pl ^ Rqa {u)g>^ Rqla {uj) 
VRVRgt,R q MGd<,,wM + VrV^^mo^M 



-v r v* r 



Q< 

Rq' a,Rqcr 



where summations over momentum indices are assumed. r.h.s. of Eq. (|C6I) can be expressed in terms of the dot 
Here, the Keldysh Green's functions which appear on the Green's functions, Q^ a da (uj) and 0^^ a (u)- 
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